Surface solitons in left-handed metamaterials 
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A theory of self-induced transparency of surface TM-mode propagating along a interface sep- 
arating conventional and left-handed metamaterials is developed. A transition layer sandwiched 
between connected media is described using a model of a two-dimensional gas of quantum dots. 
Explicit analytical expressions for a surface optical soliton in the presence of single and biexciton 
transitions, depending on the magnetic permeability of the left-handed medium, are obtained with 
realistic parameters which can be reached in current experiments. It is shown that the sign of the 
total energy flow the surface mode depends on the material parameters of the quantum dots and 
O |' the connected media. 
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The electromagnetic wave propagation through media is characterized by the electric permittivity s, and the 
magnetic permeability /z, describing the coupling of the medium to the electric and magnetic fields of the wave. The 
index of refraction of the medium n (n 2 = efi) is determined by the dispersion law for the electromagnetic plane 
wave. For plane waves propagating in isotropic regular media having simultaneously positive quantities e and /x, the 
vectors of the electric E and magnetic H fields and wave vector k constitute a right-handed triplet of vectors. In such 
media the Poynting vector S, which is the cross product of the vector of the electric and the magnetic fields, does 
coincide with the wave vector. Such media are sometimes labeled as the ordinary right-handed materials (RHM). 
In conventional media the quantities e and /i are positive. But Veselago predicted that plane electromagnetic wave 
propagation in an isotropic medium having simultaneously negative permittivity e < and permeability /i < have 
several other dramatically unusual electrodynamic properties compared with regular RHM [1]. While the Poynting 
vector for a plane wave still gives the direction of energy flow, the wave vector will point in the opposite direction 
of S and vectors E, H and k shall form a left-handed orthogonal set of vectors. In this case, they are termed 
t-H ■ left-handed materials (LHM) . Among unusual properties of the electromagnetic wave propagation in such media is 
negative refraction if the wave refracted on the same side of the surface normal as the incoming wave at an RHM/LHM 
interface. Therefore the LHM is also called a negative refractive index material (NIM). These artificial metamaterials 
\ are fascinating because they exhibit properties that do not occur in nature [2-4]. 

Pendry showed that a combination of split-ring resonators and metallic wires can lead to a negative index of 
refraction [5]. The main difficulty is to obtain a negative permeability /i < 0, which does not occur in natural 
materials. Starting with first realizations in the microwave region [6], recently, using silver as a constituent material, 
the achieved demonstrate a negative index of refraction at the red end of the visible spectrum (780 nm wavelength 
O. )[7]. 

Negative refraction could occur for surface electromagnetic waves (SEW) propagation in layered structures. Surface 
modes can propagate in a frequency region where the permittivities of two connected media have the opposite signs, 
i.e. SEW propagate at the interface between media with different " Tightness". Lately SEW have attracted much 
interest in nano-optics [8] , diverse applications [9] and also can be used as a convenient tool for studying properties of 
the structured surfaces considered as optical metamaterials [10]. Characterized peculiarity of these waves are strong 
enhancement and spatial confinement of the electromagnetic field of the wave near to the interface. In NIM interest 
to SEW is connected with the strong impact of the SEW on the image resolution of an LHM flat lens [11]. The 
properties of the SEW propagating at the interface between conventional and LHM as well as between two different 
LHM have been investigated for a intensity-dependent dielectric permittivity [12-14]. 

Another type resonant nonlinear waves can be formed within the McCall-Hahn mechanism, where a nonlinear co- 
herent interaction takes place via Rabi-oscillations of the carrier density, if the conditions for self- induced transparency 
(SIT) are fulfilled: 

wT»l, T«T h2 , (1) 

in attenuator medium. The steady-state 27r -pulse (soliton) is generated, when > ir [15]. Here T and u) are the 
width and the carrier frequency of the pulse, T\ and T2 are the longitudinal and transverse relaxation times of the 
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atomic systems, is the area of the pulse. In conventional media SIT for SEW for atomic systems and semiconductor 
quantum dots (SQD) have been investigated [16-18]. SQD, also referred to as zero-dimensional atoms(artificial atoms), 
are nanostructures which allow confinement of the charge carriers in all three spatial directions, resulting in atomic- 
like discrete energy spectra with strongly enhanced carrier lifctimes[19]. Such features make quantum dots in many 
respects similar to atoms. SIT in SQD have been investigated experimentally [20] and theoretically [21,22] for plane 
waves and waveguide modes. 

The purpose of the present Letter is to theoretically investigate the processes of the formation of surface optical 
solitons under SIT conditions in a resonance transition layer of two-dimensional sheet of inhomogeneously broadened 
SQD sandwiched at the LHM/RHM interface. 

We consider the formation of an optical soliton for surface TM-modes in a SQD transition layer sandwiched between 
conventional RHM and the LHM interface in the case when a optical pulse with width T and frequency u> » T^ 1 
is propagating along positive direction of the z axis. On the flat border of the division (x=0) between two isotropic 
media, a thin transition layer with thickness d containing a small concentration of SQD having the polarization 
P(x, z, t) = e p p(z, t)S(x), where e p is the polarization unit vector along z axis. Semi-spaces are divided at x > and 
x < 0, with conventional RHM (medium I) and LHM (medium II) with different electric permittivities £i(w) > 
and e 2 {ui) < and the magnetic permeabilities Hi(lu) > and H 2 {ui) < , respectively. In optical region of spectra 
d « A, where A is length of the surface optical wave [19]. Therefore, SIT can be modeled as SEW propagation along 
the interface LHM/RHM and infinite small thickness transition layer (monolayer) of SQDs. 

For a surface TM-mode the electric field E(E x ,0, E z ) lies in the xz plane perpendicular to the boundary of the 
division between two connected media and the magnetic field H(0, H y , 0) is directed along the axis y. The quantity 

U!(x,z,t) = [ U 1 (n,Q)e KlX e^ Qz - nt UndQ, for x<0, 



U 2 {x,z,t) = j U 2 (n,Q)e- K2X e l( - Qz - m) dndQ, for x > 0, (2) 

where 
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k? = Q 2 -Siityviifl)^-, i = 1,2. (3) 

or a Fourier-decomposition of the fields. The functions Ui :2 stand for the components (E x ,E z ,H y ) in both connected 
media. We assume translational invariance in the y-direction, so that all field quantities do not depend from the 
coordinate y, i.e. — > 0. 

Taking into account the surface current caused by presence of the SQDs, the boundary conditions for surface waves 
at x = read [16,18,23]: 

Ait dp 

H 2 , v - Hi tV = — — , E ljZ =E 2 , z . (4) 

Using the eqs.(2)-(4)we obtain the nonlinear wave equation for the E z component of the strength of electrical field at 
x = 0, in the following form[18]: 

J /(O, Q)E(Q, Q)e tiQz ~ nt UndQ = -4tt P (z, t), (5) 

where 

/(fi,Q) = £! + £!, (6) 

K 2 Kl 

E\ tZ (n, Q) — E 2tZ (Cl, Q) — E(Cl, Q). This equation is valid for any dependence of the polarization of the SQD p(z, t) 
on the strength of electrical field at x = 0. In order to determine the dependence of the polarization p(z, t) on the 
strength of electrical field at x = 0, we have to consider the structure of the energetic levels of the SQD and the 
details of the nonlinear interaction of the surface pulse with the SQD. We assume that the SQD can be described by 
the ground |0 >, exciton |2 > and biexciton |3 > states. 
The Hamiltonian of the system [18,21,22] 



H = H Q + V, 
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where 

H = Tiuj 12 \2 >< 2| + huj 13 \3 X 3| 

is the Hamiltonian of the single-exciton and biexciton states and V = —PE is the Hamiltonian of the light-quantum 
dots interaction, h is the Planck's constant, are frequencies of excitations between energetic levels i and j (i, j, = 
1,2,3). 

In general, the detunings from the resonance W13 — U12 — ui and W12 — oj, which describe the SQDs, are different. 
Under the assumption of off-resonant excitation with a constant detuning W13 — oji 2 —lu~ LO12 — L0 = A the polarization 
which is determined by intcrband transitions occurring in the quantum dots between the three energetic levels 

P = N J <?(A)[ Ml2 p 21 (A) + p 23 p 32 {A)}dA + c.c, (7) 

where TV is the uniform quantum dot density, p,i 2 = ^i 2 e, M23 = ^236*; pi 2 and ^23 are the dipole elements for the 
corresponding transitions, e is the polarization unite vector along E. We assume the dipole moments to be parallel 
to each other and to be directed along axis z, p,\ 2 — p, 23 ; s(Aw) is the inhomogeneous broadening lineshape function 
resulting from dots size fluctuations. The quantities pij are the matrix elements of the density matrix p are determined 
by the Liouville equation 

th ^r = nm > Plm ~ Pnl < mm >); 

where n, m, = 1, 2, 3. 

The solution of this equation we can find following the way presented in the work [18]. We can simplify Eq. (5) 
using the method of slowly changing profiles. For this purpose, we represent the functions E and p in the form 

E= s £ j E l Z l - p = p 12 Np Z 1 + c.c. (8) 
l=±i 

where Ei and p are the slowly varying complex amplitudes of the optical electric field and the polarization, Zi = 
e ii{kz-uit) _ rp Q g uaran ^ ee E is a real number, we set Ei = E^ t = E. This approximation is based on the 

consideration that the envelopes E vary sufficiently slowly in space and time as compared with the carrier wave 
parts-i.e., 

f)E f)E 

\ — \« u \E\, \ — \«k\E\, (9) 

and is called the slowly varying envelope approximation [15]. 

Substituting the equations (6) and (8) in the wave equation (5), and to take into account the explicit form of the 
envelope of the polarization (7) which determined from the Liouville equation, after divided the real and imaginary 
parts of the equation (5) we obtain dispersion law for surface pulse propagating on the interface between left-handed 
and right-handed media 

^2 _ ^ 2 £2^1 £2^1 ~ £ lM2 

c 2 £2 + £1 £2 - £1 

and a nonlinear wave equation in the form: 

(f)= = T-^-^. (11, 
where the width of the pulse T is determined by the equation 



T = TTT^tW / ff(A')F(A')dA', (12) 



c is the speed of light in vacuum, ( = t — y , V is the constant pulse velocity. 

/q ,/ _ df £i{uj) £2(^)1 



v„ = - 



jr, fh - QQ\n^,Q=k - fc [^ + ^3" 
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^-^b= W ,Q= fe -^-{^ + ^ ^2 [2eiW+fi-^-]}|n=u,,Q=fc, (13) 



where 



w 2 

K 2 (w, fc) = fc 2 - £i{u))lli{u)-^ , (14) 

w and fc are frequency and wave number of the carrier wave, v g is group velocity of linear SEW. 
The solution of Eq.(ll) for the envelope function has the form [15,18] 

E=^ f sech t -^±, (15) 

where = and we can determine the constant velocity of the 2tt pulse of the SEW: 

I l_ AnNfa f g(A')dA' 

V v g + Kf>,v g J T- 2 + A' 2 - 1 j 

In optical frequency region the magnetic permeability loses its usual physical meaning, in RHM we must put 
magnetic permeability /Lti = 1, otherwise would be an over- refinement [23]. We have to note that the description of 
an isotropic NIM in terms e(w) and fi(w) is not unique. Besides the e(u>) and approach for NIM, there is also 
an alternative description, based on the generalized, spatial dispersive permittivity i(u,k). In this approach, the 
non-local function e(cu,k), depending besides u> also on the wave vector fc, describes both electrical and magnetic 
responses of the medium and different optical effects in NIM [24] . 

The Eqs. (2),(10),(12)-(16) determine the parameters of the surface soliton for any value of x,z and t and show 
that for the existence of a soliton it is necessary that conditions f^>0 and V < v g arc fulfilled. Parameters of the 
surface optical solitons depend not only on the SQD parameters and the permittivites of the two interface media, but 
(unlike the case of interface of the two conventional connected media [18] ) also depends on the magnetic permeability 
ju 2 (w) and its derivative ^f^|n=u, of the LHM. 

The time-averaged over a period 2tt/uj of oscillation of the field Poynting vector of the TM mode (S), which is 
associated with the energy flow of the pulse, has the z components in the RHM-medium 1 and LHM- medium 2 in 
the following form[13, 25,26] 

<*-> - ^k^ v ""' - i^i & i v ~- (17 » 

The corresponding total energy flow 

N = f (S hz )dx + [ +OC (S 2 , z )dx = 4^(— + — )" ( 18 ) 

J-oo Jo 167TW eiKl £2K 2 

To take into account that in LHM e 2 — — 1£~2 1 the quantities (Si z ) and (S2 tZ ) have opposite signs and consequently 
the surface TM-mode resonance soliton will be have a vortexlike distribution of the Poynting vector [13,25]. 

Equation (18) has the general form and valid for any boundary condition for the strength of magnetic field of the 
SEW. The total energy flow is positive and coincide with the wave vector when the condition is satisfied: 

ki £i \H 2 \ 2 < K 2 \e 2 \\Hi\ 2 . (19) 

In case when the transition layer is absent, the boundary condition for envelopes of the strength of magnetic field of 
the surface TM-mode has the following form: H 2 = Hi, and because for SEW si < \e 2 \ and Ki < k 2 , the condition 
(19) always is satisfied and total energy flow N is always positive [25,26]. 

The situation is changed when the resonance transition layer is included. From the boundary condition (4) the 
connection between H 2 and Hi has following form 

\H 2 \ 2 = \Hi\ 2 + tR(p*Hi-pH*)+R 2 \p\ 2 (20) 

where 

Attuj^i 2 N 
H — . 
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From equation (20) it is evident that the resonance transition layer cause the change the total energy flow of the 
surface TM-mode and it will be positive or negative depending from the parameter R and polarization p. When the 
condition 

KiS^H^ + iRir^-pH^ + R 2 ^ 2 } < K 2 | £2 ||iii| 2 (21) 

is fulfilled the total energy flow will be positive (forward SEW), otherwise negative (backward SEW). 

In the considered case, the quantities H 2 and H\ are real functions and hence the equations (20) is simplified 

me^H 2 +iR(p* -p)H 1 +R 2 \p\ 2 } < K2 \e 2 \H 2 . (22) 

Consequently, when the resonance transition layer is included, depending from the parameters of the SEW, SQD 
and the connected LHM/RHM, the total energy flow can have a positive or negative sign as well. This property is 
general and are valid for any transition layer which causes the surface polarization (current) and consequently change 
the boundary conditions. This statement are valid for any linear or nonlinear wave processes for TM-mode in NIM 
and also for SEW propagating along interface between conventional media. Under the condition TV = we obtain the 
localization conditions of the surface excitations. 

In the optical NIM precious metallic nanostructures are used, therefore the losses are significant for the surface 
optical waves in LHM. But in present work the questions regarding of the losses of the surface optical modes and 
their amplifications are not considered (see, for example [27,28]). In metals and LHM the SEW occurs when the 
carrier frequency u> is below of plasma frequency u> p . For the coherent interaction of the pulse of surface wave with 
medium the duration of the pulse should be much shorter than the characteristic plasmonic oscillation damping time 
2ir/"f. Silver is known as significantly lower losses than other metals at the optical frequencies. The plasma frequency 
for silver lo p = 2ir x 2, 18 • 10 15 s _1 and damping rate is 7 = 2tt x 5, 08 • 10 12 [3]. On the other hand a planar 
array of split-ring resonators can be fabricated upon gallium arsenide (GaAs) substrate [29]. On the surface or at 
the boundary of semiconductors (for instant GaAs or InAs) with another medium, a small concentration of SQD can 
be grown. Soliton of the SEW in regular media when one of the connected medium was semiconductor with SQD is 
investigated [18]. SQD is one of the promising object also for amplification of the optical waves [28,30]. Transverse 
relaxation times of the quantum dots, which are of order of nanoseconds to several tens of picoseconds [31,32], is longer 
than the characteristic plasmonic oscillation damping time 2-71-/7. The envelope approximation (9) is appropriate for 
pulses when the spectral width of the pulse much smaller than the u and valid for pulses with width T > 20/s [33] . 
Consequently pulse with such width which acceptable for investigation of the nonresonance nonlinear waves in the 
LHM [12,13], is satisfying also the eq.(l) the conditions of SIT in SQD . This circumstances allow to hope that SIT 
in LHM/RHM with SQD transition layer can be experimentally observable. 

In conclusion, we considered case when a transition layer contains SQD, however, the presented results have a 
more general character and are valid also in case when the transition layer contains usual optical impurity atoms 
sandwiched between LHM/RHM. 
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